Inferring introduction routes of invasive species using approximate Bayesian computation on microsatellite data

Determining the routes of introduction provides not only information about the history of an invasion process, but also information about the origin and construction of the genetic composition of the invading population. It remains difficult, however, to infer introduction routes from molecular data because of a lack of appropriate methods. We evaluate here the use of an approximate Bayesian computation (ABC) method for estimating the probabilities of introduction routes of invasive populations based on microsatellite data. We considered the crucial case of a single source population from which two invasive populations originated either serially from a single introduction event or from two independent introduction events. Using simulated datasets, we found that the method gave correct inferences and was robust to many erroneous beliefs. The method was also more efficient than traditional methods based on raw values of statistics such as assignment likelihood or pairwise F(ST). We illustrate some of the features of our ABC method, using real microsatellite datasets obtained for invasive populations of the western corn rootworm, Diabrotica virgifera virgifera. Most computations were performed with the DIYABC program (http://www1.montpellier.inra.fr/CBGP/diyabc/)

Anarchy in the UK: Detailed genetic analysis of worker reproduction in a naturally occurring British anarchistic honeybee, Apis mellifera, colony using DNA microsatellites

Anarchistic behaviour is a very rare phenotype of honeybee colonies. In an anarchistic colony, many workers’ sons are reared in the presence of the queen. Anarchy has previously been described in only two Australian colonies. Here we report on a first detailed genetic analysis of a British anarchistic colony. Male pupae were present in great abundance above the queen excluder, which was clearly indicative of extensive worker reproduction and is the hallmark of anarchy. Seventeen microsatellite loci were used to analyse these male pupae, allowing us to address whether all the males were indeed workers’ sons, and how many worker patrilines and individual workers produced them. In the sample, 95 of 96 of the males were definitely workers’ sons. Given that ≈ 1% of workers’ sons were genetically indistinguishable from queen’s sons, this suggests that workers do not move any queen-laid eggs between the part of the colony where the queen is present to the area above the queen excluder which the queen cannot enter. The colony had 16 patrilines, with an effective number of patrilines of 9.85. The 75 males that could be assigned with certainty to a patriline came from 7 patrilines, with an effective number of 4.21. They were the offspring of at least 19 workers. This is in contrast to the two previously studied Australian naturally occurring anarchist colonies, in which most of the workers’ sons were offspring of one patriline. The high number of patrilines producing males leads to a low mean relatedness between laying workers and males of the colony. We discuss the importance of studying such colonies in the understanding of worker policing and its evolution

A Yule-Simon process with memory

The Yule-Simon model has been used as a tool to describe the growth of diverse systems, acquiring a paradigmatic character in many fields of research. Here we study a modified Yule-Simon model that takes into account the full history of the system by means of an hyperbolic memory kernel. We show how the memory kernel changes the properties of preferential attachment and provide an approximate analytical solution for the frequency distribution density as well as for the frequency-rank distribution.Comment: 7 pages, 5 figures; accepted for publication in Europhysics Letter

Individual phenotypic variation reduces interaction strengths in a consumer–resource system

Natural populations often show variation in traits that can affect the strength of interspecific interactions. Interaction strengths in turn influence the fate of pairwise interacting populations and the stability of food webs. Understanding the mechanisms relating individual phenotypic variation to interaction strengths is thus central to assess how trait variation affects population and community dynamics. We incorporated nonheritable variation in attack rates and handling times into a classical consumer–resource model to investigate how variation may alter interaction strengths, population dynamics, species persistence, and invasiveness. We found that individual variation influences species persistence through its effect on interaction strengths. In many scenarios, interaction strengths decrease with variation, which in turn affects species coexistence and stability. Because environmental change alters the direction and strength of selection acting upon phenotypic traits, our results have implications for species coexistence in a context of habitat fragmentation, climate change, and the arrival of exotic species to native ecosystems

Chloroplast microsatellites: measures of genetic diversity and the effect of homoplasy

Chloroplast microsatellites have been widely used in population genetic studies of conifers in recent years. However, their haplotype configurations suggest that they could have high levels of homoplasy, thus limiting the power of these molecular markers. A coalescent-based computer simulation was used to explore the influence of homoplasy on measures of genetic diversity based on chloroplast microsatellites. The conditions of the simulation were defined to fit isolated populations originating from the colonization of one single haplotype into an area left available after a glacial retreat. Simulated data were compared with empirical data available from the literature for a species of Pinus that has expanded north after the Last Glacial Maximum. In the evaluation of genetic diversity, homoplasy was found to have little influence on Nei's unbiased haplotype diversity (H(E)) while Goldstein's genetic distance estimates (D2sh) were much more affected. The effect of the number of chloroplast microsatellite loci for evaluation of genetic diversity is also discussed

Genomic footprints of a biological invasion: Introduction from Asia and dispersal in Europe of the topmouth gudgeon (Pseudorasbora parva)

Facilitated by the intensification of global trading, the introduction and dispersal of species to areas in which they are historically non-native is nowadays common. From an evolutionary standpoint, invasions are paradoxical: not only non-native environments could be different from native ones for which introduced individuals would be ill-adapted, but also small founding population size should be associated with reduced adaptive potential. As such, biological invasions are considered valuable real-time evolutionary experiments. Here, we investigated the population structure and adaptive potential of the highly invasive topmouth gudgeon (Pseudorasbora parva) across Europe and East Asia. We RAD-sequenced 301 specimens from sixteen populations and three distinct within-catchment invaded regions as well as two locations in the native range. With 13,785 single nucleotide polymorphisms, we provide conclusive evidence for a genome-wide signature of two distinct invasion events, in Slovakia and Turkey, each originating from a specific area in the native range. A third invaded area, in France, appears to be the result of dispersal within the invasive range. Few loci showed signs of selection, the vast majority of which being identified in the Slovakian region. Functional annotation suggests that faster early stage development, resistance to pollution and immunocompetence contribute to the invasion success of the local habitats. By showing that populations in the invasive range have different evolutionary histories, our study reinforces the idea that populations, rather than species, are the units to consider in invasion biology

Secondary contact and admixture between independently invading populations of the Western corn rootworm, diabrotica virgifera virgifera in Europe

The western corn rootworm, Diabrotica virgifera virgifera (Coleoptera: Chrysomelidae), is one of the most destructive pests of corn in North America and is currently invading Europe. The two major invasive outbreaks of rootworm in Europe have occurred, in North-West Italy and in Central and South-Eastern Europe. These two outbreaks originated from independent introductions from North America. Secondary contact probably occurred in North Italy between these two outbreaks, in 2008. We used 13 microsatellite markers to conduct a population genetics study, to demonstrate that this geographic contact resulted in a zone of admixture in the Italian region of Veneto. We show that i) genetic variation is greater in the contact zone than in the parental outbreaks; ii) several signs of admixture were detected in some Venetian samples, in a Bayesian analysis of the population structure and in an approximate Bayesian computation analysis of historical scenarios and, finally, iii) allelic frequency clines were observed at microsatellite loci. The contact between the invasive outbreaks in North-West Italy and Central and South-Eastern Europe resulted in a zone of admixture, with particular characteristics. The evolutionary implications of the existence of a zone of admixture in Northern Italy and their possible impact on the invasion success of the western corn rootworm are discussed